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Abstract: A review of results from lattice studies using improved and scale 
controlled cooling methods is presented and their significance is discussed. 
The improvement of the action ensures stable instanton solutions of physical 
sizes. The scale controlled cooling can be generally used as a gauge invariant 
low pass filter to extract the physics from noisy MC configurations; in 
particular it preserves instanton-antiinstanton pairs selected according to 
their interaction. We apply these methods to analyze various features of 
the topological structure of the Yang-Mills vacuum in a scaling invariant 
way. 

1. Introduction 

This article presents results from a series of works involving besides the 

above author the following persons in various combinations: 

PH. DE FORCRAND (ETH Zurich and TH-Div., CERN, Switzerland), 

M. GARCIA PEREZ (Dept. Ffs. Teor., Univ. Aut. de Madrid, Spain), 

T. HASHIMOTO (Dept. of Appl. Physics, Fukui Univ., Japan), 

JE. HETRICK (Physics Dept., U. of the Pacific, Stockton, USA), 

S. HIOKI (Dept. of Physics, Tezukayama Univ., Nara, Japan), 

E. LAERMANN (FB Physik, Univ. Bielefeld, Germany), 

JF. LAGAE (HEP Div., Argonne Nat. Lab., USA), 

H. MATSUFURU (RCNP, Osaka, Japan), 

O. MIYAMURA, T. UMEDA (Dept. of Physics, Hiroshima Univ., Japan), 

A. NAKAMURA (RUSE, Hiroshima Univ., Japan), 

O. PHILIPSEN (Th. Div., CERN, Geneve, Switzerland), 

T. TAKAISHI (Hiroshima University of Economics, Japan). 

The simulations and results on which the following discussions are based 

are described in [1], [2], [3], [4], [5], [6], [7]. 
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2. Physical Questions and Problems of the Analysis 

Instantons are self-dual solutions of the Euclidean field equations map- 
ping the SU{2) color group space on the 3-d sphere at spatial infinity (for 
SU(N C ) the SU(2) subgroups are relevant). They are defined among other 
by winding number (topological charge) and position and size parameters. 
The R 4 (anti-)instanton (A)I has charge (q) and action (s) densities [8]: 



q[X) 8tt 2 * ttV 4 




, s(x) = S \q(x)\, (1) 



Q = J d 4 xq(x) = iinteger, S = J d 4 xs(x) = |Q|So, 5*0 



8tt 2 



In T 4 only \Q\ ^ 1 exist [9]. In R 3 T (relevant at non-zero temperature T) 
the "caloron" [10] describes objects flattening in the Euclidean time direc- 
tion and of limited spatial size p < \^3/irT. Finally, various kinds of twisted 
boundary conditions allow for new solutions [11], [12]. Superpositions of N 
I's or A's also correspond to (higher) minima of the action: S = NSq. A 
pair I-A is not, however, a minimum and S IA depends also on "overlap" u>: 



S lA = 2S - < 2S , u}=(jn + PA)/d 1A , djA = |x? - x° A \. (2) 

Effects associated with instantons are (for a review see [13]): 

- the Ua(1) symmetry breaking (via the Witten - Veneziano formula[14]), 

- chiral symmetry breaking (via zero modes of the Dirac operator) , 

- dynamical effects for the physics at intermediate distances. 

The first one only involves global topological properties (susceptibility). 
The other possible effects depend on details of the local structure, like den- 
sity and size distribution of instantons, and the calculation involves models. 
In the instanton liquid model [15], for instance, which is the standard model 
of instanton physics, various predictions have been made [13]. Numerical 
simulations should test the ingredients used in these models. Thereby: 

- The instantons can be in a gas, liquid, or crystalline phase. The diluteness 
is expressed by the "packing fraction" 

/ = (3) 
J 2V 

- The size distribution controls diluteness and I-R properties. There is no 
prediction for large sizes. For small sizes the dilute gas approximation gives: 



P{p)~P p , dilute = -5 + b, 6=llAT c /3 



(4) 
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- Denoting by Nj(Na) the number of I's (A's) in a configuration, and with 
N = Ni + Na, Q = Ni - N A we write: 

cp =(<iV 2 >- (Nf) /(N). (5) 

For a poissonian distribution cp = 1, while low energy sum rules suggest 
cp = | [16]. A further measure for the character of the I, A distribution is: 

r WVa) - (Ni)(N A ) _ (iV 2 ) - (Q 2 ) - (AQ 2 

mt ^((n?) - (N imN i) - { N A )i) m + m-(N)i u 

which tests the correlation between I's and A's: c- ul t is if there is no such 
correlation, but goes to 1 if I and A appear correlated (e.g., pairwise). 

For non-perturbative analysis one resorts to lattice simulations. The 
Monte Carlo method is sometimes considered to be a tool which produces 
numerical results, without improving our understanding of the physical 
mechanisms. This view, however, has turned out to be too narrow and dif- 
ferent structural questions have been investigated by means of Monte Carlo 
calculations. The problems such an analysis of the topological structure en- 
counters are twofold. Problems of principle are: 

- UV-divergences in correlation functions like (q(x)q(0)) which contain 
the susceptibility information in a contact term (notice that (q(x)q(0)) < 
for disjoint supports, e.g. for x > 2a if F^ u is defined on plaquettes) [17]. 

- Dislocations: concentration of topological charge on one or few pla- 
quettes. They have action smaller than Sq and may be produced copiously, 
spoiling the measurement of Q and the continuum limit [18]. 

- Since Sj^j. < 2Sq close pairs can be easily produced. With increasing 
overlap it becomes difficult to define genuine I-A pairs as distinguished from 
short range density fluctuations and the description in terms of I's and A's 
becomes questionable [3]. 

It appears therefore that a certain amount of smoothing may be neces- 
sary to define topology on the lattice. Practical problems are: 

- The need to develop smoothing methods which can be controlled in 
a physical way. They should be based on independent physical arguments, 
leave the investigated structure undistorted and allow to observe its prop- 
erties without making assumptions about them. In particular scaling be- 
haviour of topological distributions should not be assumed but revealed. 

- The description of the topological structure in terms of I's and A's. 
Overlapping structures, especially with opposite charges, become difficult 
to fit with a superposition ansatz based on (1). One needs to introduce 
various criteria and some measure for the adequacy of the description which, 
however, remains imperfect. 
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3. Scale Controlled Cooling [5], [6] 

The problem of developing good smoothing methods has been approached 
in a number of ways: standard cooling, smearing, blocking, combination of 
them and others [19], [20], [21], [22]. The smoothing is controlled by monitor- 
ing, by giving a blocking step or by tuning smearing to approximate a RG 
transformation with fixed point action etc. The most natural control for 
a smoothing method, however, is by giving the smoothing radius or wave 
length threshold and we should like, therefore, to develop a method sys- 
tematically controlled by a physical scale. We start from cooling, since this 
procedure is analytically defined (by reference to the equations of motion) 
and has as fixed point nontrivial semi-classically relevant configurations 
with N = \Q\. Cooling, which proceeds by local minimization of a given 
lattice action, works as a diffusion process, smoothing out increasingly large 
regions [23]. Thereby properties of the physical spectrum measured from 
correlations at a given distance are affected, in particular the string ten- 
sion drops rapidly with the number of cooling sweeps n c [24]. Instantons 
may be distorted and lost during cooling by bad scaling properties of the 
action and also by I- A annihilation. The susceptibility may be spoiled by 
dislocations. We thus need: 

1) To use an action with practically scale invariant instanton solutions and 
a sharp dislocation threshold fixed in lattice units to eliminate UV noise 
(dislocations) while shrinking in approaching continuum, and 

2) To ensure that cooling is controlled by a physical smoothing scale such 
that the structure at larger scales remains unchanged. Monitoring and en- 
gineering should not be necessary, since this may introduce arbitrariness. 

The Restricted Improved Cooling (RIC) introduced in [5] fulfills these re- 
quirements. RIC acts as a gauge invariant low pass filter preserving physics 
at scales above a cooling radius r c which can be fixed unequivocally, while 
smoothing out the structure below r c . In particular the string tension is 
preserved beyond r c , instantons are stable, dislocations are eliminated, and 
I-A pairs are retained above an overlap threshold depending on r c . 

RIC uses the action of the Improved Cooling algorithm (IC) with 5 pla- 
nar, fundamental Wilson loops [25], [1]. This action is correct to order C(a 6 ) 
and has a dislocation threshold po/a = po ~ 2.3, below which short range 
topological structure is smoothed out (note that po — > in approaching con- 
tinuum). Above po, instantons are stable under IC to any degree of cooling 
(however, I-A pairs annihilate). The corresponding improved charge density 
using the same combination of loops leads to an integer charge already after 
a few cooling sweeps and stable thereafter. See Fig. 1. 

Recall that the cooling algorithm is derived from the equations of motion 

V^W^x? - W^{x)U^ = 0, (7) 
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Figure 1. IC and RIC properties: Evolution of sizes of single instantons (upper left). 
Improved (5Li) and Wilson action vs instanton size given in lattice units; p pea k is defined 
in [1] (upper right). Evolution of the improved charge Q for 20 5(7(3) MC configurations 
[7] (lower left). Saturation of cooling for RIC: fraction of updated links as function of the 
cooling sweeps - for stopped standard cooling or IC this would be a step function (lower 
right). 

where W is the sum of staples connected to the link U^(x) in the action, 
as (we restrict here to SU(2), for the general case see [5]): 



U U' = V = W/\\W\\, \\W\\ 2 = -Ty(WW ] ). 



(8) 



We define RIC by the constraint that only those links be updated, which 
violate the equation of motion by more than some chosen threshold 1 : 



U V iff A M (x) 2 = a- 6 Tr(l - C/W) > 5 2 . 
1 We thank F. Niedermayer for this suggestion. 



(9) 



6 



We have A^(x) oc —Tr((D L ,F Ufl (x)) 2 ) in continuum limit [26]. Thus 5 con- 
trols the energy of the fluctuations around classical solutions and acts as 
a filter for short wavelengths. Since it uses the same action RIC has the 
same scaling properties as IC. However, since RIC does not update links 
already close to a solution, it changes fewer links after every iteration until 
the algorithm saturates - see Fig. 1. (9) defines a constrained minimization 
and the smoothing is homogeneous over the lattice (see [5]). 

Since the parameter 6 which defines the cooling is already a physical 
quantity it should be related to a physical scale. For Yang Mills theory the 
latter involves the string tension a. We calculate the "effective mass" M(t) 
from correlation functions of spatial Polyakov loops separated by t steps in 
time. Asymptotically M(t) ~ N s a up to finite size corrections. 




Figure 2. Left and middle: M(t) for the a = 0.12 (0.06) fm lattice [5]. Horizontal bands 
represent standard results for a. Right: The smoothing scale of RIC r c (S) vs <5~ 1//3 for 
the a = 0.12 fm (diamonds) and a = 0.06 fm (squares) lattices. Also shown is the fit (10). 

We present in Fig. 2 M(t) for SU(2). Defining r c {5) as the distance t 
at which M(t) on (5-cooled configurations starts to agree with the uncooled 
value (obtained by fuzzing and fitted to a smooth function of t) we arrive 
at the results in Fig. 2 (right) showing the scaling behaviour expected from 
the continuum limit of the parameter 5, compatible with 

r c (5) ~ 0.8 <T 1/3 . (10) 

The behaviour of instantons is similar under RIC and IC. RIC has a 
dislocation threshold depending on a 3 5, po(a 3 5) = po(a' i 5)a approaching 
the IC dislocation threshold from below: po(a 3 5) < po(0) = po = 2.3. Fix- 
ing S and rescaling (3 to approach the continuum limit the RIC threshold 
shrinks similarly to the IC one. Instantons above the dislocation threshold 
are preserved unchanged both by RIC (independently on r c ) and by IC (in- 
dependently on cooling sweeps), since both algorithms use the same action 
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with scale invariant instanton solutions. I-A pairs, however, do not annihi- 
late under RIC the way they do for the other cooling methods (including 
IC): since the distortion of the partners in a pair depends on the overlap, 
RIC preserves pairs below some overlap-threshold. More precisely, there is 
a well defined relation between A M (measured at the center of the partners 
or at mid-distance between them) and Sf^ eq. (2). Thus RIC stabilizes I-A 
pairs with Sj^j. below a threshold depending on r c (<5). 

4. Topological Properties of Yang-Mills Theories at T = 

4.1. GLOBAL PROPERTIES [1],[3],[5],[6] 

The topological charge given by IC is already after a few cooling sweeps 
an integer to about 1% and stable thereafter. Similarly, Q given by RIC is 
also an integer and independent on S if the cooling radius is larger than the 
dislocation threshold, r c > po ~ 2a. The topological charge distribution is 
therefore well defined. It also agrees excellently, e.g., with the distribution 
obtained by using overlap fermions [27]. The topological susceptibility ob- 
tained in this way is well defined, shows good scaling properties and agrees 
with the Witten-Veneziano relation. See Table 1. 
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0.115 


0.103 


0.0855 


X 1/4 (MeV) 


195(4) 


200(8) 


184(6) 


182(7) 


134(10) 


102(5) 






TABLE 1. x 1/4 for pure Yang-Mills at T = [1],[3],[5] (average) and for QCD at 
T ~ 0.85, 0.98 and 1.13T C [3]. The results are given for n c > 20 or r c > 2a. 



4.2. LOCAL PROPERTIES [1],[3],[5],[6],[7] 

4.2.1. Description of the I-A Ensemble and Overview of the Results 
In extracting the instanton information out of the configurations we ap- 
proximate the action and charge density by a superposition of self-dual, 
non-interacting I's and A's parameterized through the 1-instanton BPST 
ansatz (1) (partially corrected for periodicity by adding images). We mea- 
sure the departure of the real action and charge density from the above 
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non-interacting ansatz through the quantities 




Jd 4 x\s(x) - s fit (V)| 2 
/ d 4 x\s(x)\ 2 




Jd 4 x\q(x) - gfit(x)| 2 
Jd 4 x\q(x)\ 2 



(11) 



In the actual proceeding every time an instanton candidate is located, it is 
counted only if, by adding it to sgt and gg t , e S)9 simultaneously decrease. 
Of course this description misses more complicated topological structures 
(e.g., fractional charge) - it is a question of self-consistency to apply it 
within a picture of the instanton ensemble consisting of typical, self-dual 
I's and A's in weakly interacting superpositions. Attempts to use a more 
general description have been made in the literature [28], their theoretical 
understanding appears to us difficult, however. 

This analysis is performed at a given degree of smoothness and we ask 
about the stability of the results with the latter. Our attitude is that signif- 
icant dependency on the smoothing degree can reflect, may be unwanted, 
but real properties of the instanton ensemble. To judge about this scaling 
tests are essential. Normal cooling and even IC reveal the smoothing depen- 
dence as dependence on the number of cooling steps. The latter, however, 
has only an indirect physical meaning [23]. It needs to be calibrated in some 
way with (5 (see [1] , [31] and for a recent systematic attempt [32]), the cri- 
terion for the calibration remaining, however, just the property which one 
wants to prove, namely scaling. RIC, on the other hand, has an automatic, 
natural calibration via the physical (dimensionfull) cooling scale r c . The 
cooling is independent on (3 if r c is larger than the shrinking dislocation 
threshold po, and therefore provides a genuine scaling check [5]. We then 
consider any r c dependence which subsists the scaling test as potentially 
relevant and do not attempt a further treatment of the data (e.g., extrap- 
olation to smoothing degree [22], [29]) since we are not clear about the 
physical basis of the dependence and of the treatment. 

An overview of the T = RIC and IC SU{2) results is given in Table 
2, for SU(3) and QCD we refer to [3]. 

4.2.2. Instanton Size Distributions 

As can be seen from Fig. 3 the size distributions show remarkable stability 
under changing the RIC smoothing scale and under rescaling of the lattice 
spacing, up to a limited variation of the position of the maximum (notice 
that these are normalized distributions, the number of instantons varies 
by more than a factor 5 between small and large r c ). Also the IC data 
show the same stability, both for SU(2) and for SU(3) - see Fig. 3. The 
distributions are little (at larger a) or not affected (at smaller a) by the 
dislocation threshold. The average sizes p av are 0.4 i 0.04 fm for ( S'C/(2) and 
0.55±0.06 fm for SU(3), with full widths of about 0.2 (0.25) fm, respectively. 
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8 
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(\Q\) 
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{diA) 


e q (e s ) 
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15 
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4.21(2) 


2.91(6) 


0.54 
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0.45(0.35) 
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15 


2.89 


0.60(5) 


2.70(3) 


2.85(8) 


0.59 


0.58 


0.35(0.30) 


0.44 


15 


IC(20) 




2.06(2) 


2.91(8) 


0.61 


0.61 


0.30(0.30) 


0.42 




0.06 


23.15 


0.27(3) 


18.7(1) 


1.7(1) 


0.40 


0.30 


0.60(0.45) 


0.35 


15 


11.57 


0.33(3) 


9.6(1) 


1.7(1) 


0.48 


0.39 


0.45(0.40) 


0.38 


15 


4.92 


0.50(3) 


4.35(6) 


1.8(1) 


0.57 


0.53 


0.25(0.25) 


0.40 


15 


IC(50) 




2.96(5) 


1.7(1) 


0.59 


0.57 


0.25(0.25) 


0.37 



TABLE 2. SU(2) results for IC at 20 and 50 cooling sweeps and for RIC at various 
8 [5] - see also Table 1. The lattice spacing (a), the cooling radius (r c (<5)), distances 
and sizes are given in fm. The density (N/V) is in fm~ 4 and 8 in fm -3 . The errors 
are statistical only, when they are not explicitly quoted they amount to about 1 in 
the last indicated digit. 



Since the variation of p av , which may have different origins [30], is rather 
limited we are tempted to consider this quantity as well defined (inside 
errors accounting for the various uncertainties) - see also Fig. 5. The small 
size branch of the distribution grows more slowly than the power law of 
the dense gas approximation (a fit gives p = 1.2 ± 0.5 for SU(2) and p = 
2.5 ± 1.0 for SU(3), instead of ^dilute = 2.33 and Pdiiutc = 6, respectively) 
which indicates that in the region of small but relevant sizes the dilute gas 
picture is not adequate. The large size behaviour of the distributions cannot 
be determined from these simulations since instantons of size significantly 
above 1/3 of the lattice can show finite lattice size effects and already p av is 
near this point (0.48fm for SU{2) and 0.53 fm for S77(3); but the existence 
and position of the maximum are safe from this point of view). 

4.2.3. Instanton Density and Pair Properties 

The number N of I's and A's observed at small cooling radius is essen- 
tially due to pairs (|Q| is constant and relatively small) - see Table 2. We 
notice that while N increases by roughly a factor 5 when r c is decreased 
from about 0.6 fm to about 0.25 fm, p av only varies by about 10%. Con- 
sequently, the packing fraction / (3) also increases approaching 1, which 
indicates that at small r c we observe a rather dense ensemble, strongly 
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Left: Normalized size distributions for 5(7(2) from RIC and IC on the 
12 3 36, pbc, a = 0.12 fm lattice (squares, triangles) and on the 24 4 , tbc, a = 0.06 fm 
lattice (crosses) [5]. The RIC data correspond to r c ~ 0.56, 0.47, 0.35 and 0.28 fm 
(S = 2.89, 4.92, 11.57 and 23.15 fm~ 3 , respectively). Vertical dotted (dashed) lines in- 
dicate the IC dislocation threshold p = 2.3a for a = 0.12 (0.06) fm. The corresponding 
RIC dislocation thresholds are lower. Right: Normalized size distributions for SU(3) from 
IC on the 12 4 , tbc, a = 0.134 fm lattice (squares, circles, triangles) and on the 16 4 , tbc, 
a = 0.1 fm lattice (octagons) [3]. The curves are Gaussian fits, the dashed lines po = 2.3a. 




Figure 4- Non-normalized size and overlap distributions from RIC at a = 0.12 fm [5] 
for smoothing scales r c ~ 0.56 (squres), ~ 0.47 (pluses) and ~ 0.28 (diamonds) fm. The 
scale for sizes is 0.06 fm. 

deviating from Poisson distribution (cp eq. (5) is about 0.5) and with a 
tendency to charge neutralization in small regions (see [1],[31]). The over- 
lap (2) distribution shifts toward larger values [6] while the average average 
distance of an I to the nearest A {(Ija) decreases. See Fig. 4. The situation 
is succinctly described in Fig. 5. Here we plot SU{2) RIC data for averages 
from two lattices (a = 0.12 fm and a = 0.06 fm) at given r c = 0.85 -1 / 3 , 
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Figure 5. Smoothing scale dependence of: x 1 ^ 4 (upper left), average 1(A) size (upper 
right), average I-A distance, (lower left) and average density (lower right) for SU(2) at 
a = 0.12 fm (crosses) and a — 0.06 fm (squares) [5]. r c is taken from the fit (10), the 
horizontal errors (not shown) are about 0.04 fm - see also Table 2. The lines illustrate 
linear (dashed), r^ 3 (dashed-dotted) or r~ 4 (dotted) behaviour. 

see eq. (10), i.e., the dimensionless 8a~ 3 is simply rescaled with a at fixed 
5. As already noticed, we see that the susceptibility x an d the average size 
p av scale well with the cut-off and show no (%) or little dependence (p av ) 
on the smoothing scale. Also the average I-A distance and the density of 
1(A) 's have comparable values on the two lattices and seem therefore to 
scale. However, their dependence on r c is disturbing: the instanton den- 
sity increases as an inverse power of r c and the average I-A distance is 
practically given by r c itself. Simultaneously with increasing density also 
the quality of the fit decreases (e q , e s increase, see Table 2). This suggests 
that the structure revealed at small smoothing scale (large frequencies) is 
strongly overlapping and distorted. Generally it appears that a "typical" 
density or I-A distance could only be defined by first fixing the smoothing 
scale. By contrast, the average size seems to have a well defined meaning 
for itself. Finally, restricting smoothing to scales below 0.25 fm retains too 
much short range fluctuations. We stress that the cut-off rescaling of these 
results is done without any ad-hoc fit of cooling parameters. 2 

2 The small, systematic difference of a few percents between the "a = 0.12 fm" and the 
"a = 0.06 fm" data for all observables (including a and x) [5] may be due to a minimal 
uncertainty in assessing the value of a [1] as taken from the literature. 
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5. Topology at T > and the Chiral Connection [2] , [3] , [4] , [7] 

A first impression on the connection between topology and chiral prop- 
erties (see [33] for a recent review) is obtained at T = 0, e.g., by using 
the definition of the topological charge via the Ward identities for Wilson 
fermions [34] or via the overlap formalism [27] and comparing with the 
improved topological charge measured by IC or RIC. In Fig. 6, left plot, 
we present the topological charge distribution measured on a 12 4 , SU(2) 
lattice (a = 0.12 fm, pbc [1]) by IC and by overlap formalism [27]. See also 
[35]. 

Increasing the temperature above T c produces radical changes in the 
topological properties: the susceptibility drops and the instanton ensemble 
seems to change its character [36], [3], [2]. Above T c practically all configura- 
tions are in the Q = sector although at small r c RIC reveals a quite large 
N (a similar thing is observed with IC at small number of cooling sweeps). 
The ensemble consists of nearby opposite charges, but it is difficult to say 
if these form "molecules" or just a strongly fluctuating background. The 
description in terms of I's and A' becomes difficult (the fit deteriorates). 
On Fig. 6, right plot, we show the I-A correlation parameter c- mt (6) for the 
full QCD, MILC configurations below and above T c [37]. 

It is interesting to observe the correlation between topological and chiral 
properties in full QCD as function of the temperature. On the MILC config- 
urations [37] which use two flavours of staggered quarks we have measured 
the improved topological charge and the chirality of the lowest eigenvalues 
of the Dirac operator. Although the renormalization effects on the chirality 
are large there appears to be a strong correlation between the total chiral- 
ity carried by the small eigenvalues (unsmoothed) and the total charge of 
the configuration obtained by IC. Both are also similarly affected by the 
temperature: above T c Q = and all Dirac modes become heavy and have 
vanishing chirality. See Fig. 7 where also the I-A correlation parameter Cj nt 
(6) is drawn. 

One even observes local correlation between instantons and Dirac modes 
[2], [4], the bulks of the corresponding densities seem to overlap quite well. 

6. Conclusions 

The question: which are the vacuum excitations relevant for the dynamics 
of Yang Mills theories? can be studied in the frame of Euclidean lattice 
theory by observing the structure of the configurations generated in Monte 
Carlo simulations. Besides instantons various other proposals have been 
made, such as super-instantons [38], abelian monopoles or vortices (for re- 
cent reviews see, e.g., [39], [40], [41]). We have here been concerned with 
topological properties, that is with the problem of observing and charac- 
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Figure 6. Left: Topological charge distribution for pure SU(2) by IC at n c = 20 (di- 
amonds) and n c = 300 (crosses) and via overlap formalism (squares). The curve is an 
illustrative Gaussian with arbitrary normalization and slope corresponding to a suscep- 
tibility of (200 MeV) 4 . Right: Temperature dependence of the full QCD data of Fig. 7: 
c int for T ~ 0.85, 0.98 and 1.13T C after 20 (upper data), 50 (middle) and 150 (lower 
data) IC cooling sweeps. 
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Figure 7. QCD with dynamical quarks: Chirality /eigenvalue plots for the lowest 4 
eigenvalues using symbols corresponding to the topological charge for each configuration: 
octagon (Q — 0), cross and fancy cross (Q = ±1) and plus and fancy plus (Q = ±2). 
The three plots correspond to T ~ 0.85, 0.98 and 1.13T C [7] 

terizing the (anti-)instanton ensemble in the MC configurations. Thereby 
the I's and A's are identified as self-dual structures and fitted to the BPST 
ansatz. We obtain in this way various results (density, charge and size 
distributions etc) and we also check the scaling behaviour of these prop- 
erties under changing of the lattice spacing by factors of the order of 2. 
In agreement with other studies (see [29], [42] for recent reviews), our re- 
sults indicate that global observables like the topological susceptibility and 
charge distributions appear well defined and show good scaling properties. 
Also some local properties, like size distributions of instantons appear rea- 
sonably well defined and scaling. It seems, however, that the topological 
charge density can fluctuate strongly at small distances and a description 
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of these fluctuations in terms of instantons and anti-instantons is difficult: 
the objects do not decrease significantly in size but begin to overlap and 
become increasingly deformed. While it seems that one can give to some 
extent "physical" values (by which it is meant: invariant under rescaling 
of the coupling) for the 1(A) density and overlap if only fluctuations above 
a given wave-length are retained, these values depend on the latter and it 
appears difficult to select some "more physical" sub-ensemble characterized 
by universal parameters. 

The connection between topology and chiral properties is well observed 
at the level of the topological structure which define the topological sector: 
the instantons which appear responsible for the net topological charge ap- 
pear also strongly correlated with the Dirac modes of small eigenvalue and 
large chirality. This connection can also be followed with the temperature. 

The analysis of the MC configurations typically needs some procedure 
to get rid of short range fluctuations. We have used here a "scale controlled 
cooling algorithm", RIC, which has been developed to act as a gauge in- 
variant low pass filter. RIC smoothes out fluctuations up to a chosen scale 
given as a physical (dimensionfull) length and has no other parameter (like 
setting the number of cooling sweeps n c etc) . If the smoothing scale is taken 
large enough this algorithm reproduces standard cooling and it can be used 
under circumstances to understand the effect of the latter in its dependence 
on n c [43]. Using RIC allows for an independent scaling check, especially 
in connection with such quantities (like the I-A density) which appear to 
depend on smoothing, since the smoothing scale is independent on the lat- 
tice spacing. In general this smoothing procedure can be used to study the 
phenomena at a given physical scale and may represent therefore a good 
instrument in extracting dynamical information from MC configurations. 
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